Solar Physics 
DOI: 10.1007A 



Heliolatitude and time variations of solar wind structure 
from in situ measurements and interplanetary 
scintillation observations 

J.M. Sokol' • M. Bzowski' • M. Tokumaru^ • 
K. Fujiki^ . D.J. McComas^^ 

© Springer •••• 

Abstract The 3D structure of the solar wind and its evolution in time is needed for 
heliospheric modeling and interpretation of energetic neutral atoms observations. We 
present a model to retrieve the solar wind structure in heliolatitude and time using 
all available and complementary data sources. We determine the heliolatitude struc- 
ture of solar wind speed on a yearly time grid over the past 1.5 solar cycles based 
on remote-sensing observations of interplanetary scintillations, in situ out-of-ecliptic 
measurements from Ulysses, and in situ in-ecliptic measurements from the OMNI-2 
database. Since the in situ information on the solar wind density structure out of ecliptic 
is not available apart from the Ulysses data, we derive correlation formulae between the 
solar wind speed and density and use the information on the solar wind speed from in- 
terplanetary scintillation observations to retrieve the 3D structure of solar wind density. 
With the variations of solar wind density and speed in time and heliolatitude available 
we calculate variations in solar wind flux, dynamic pressure and charge exchange rate 
in the approximation of stationary H atoms. 

Keywords: Solar Wind: Models, Solar Wind: Observations, Radio Scintillation, Solar 
Cycle: Models, Solar Cycle: ObseiTations 

1. Introduction 

The goal of this paper is to retrieve the solar wind structure at 1 AU as a function of time 
and heliolatitude based on available in situ data sources and interplanetary scintillation 
observations for the time interval since 1990 until present that covers 1.5 solar activity 
cycles. 
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The existence of the solar wind was predicted on a theoretical basis by Parker ( 1958[ l 
and discovered experimentally by Lunnik II and Mariner 2 at the very beginning of the 
space age ( Gringauz et al. 1960t|Neugebauer and Snyder 1962| l. Regular measurements 
of its parameters began in the early 1960s and data from many spacecraft are now 
available, obtained using various techniques of observations and data processing (see 
for review |Bzowski et a/.|[2012| l. Shortly after the discovery of the solar wind (SW 
hereafter), a question of whether or not it is spherically symmetric was put forward. 
Most spacecraft with instruments to measure the SW parameters are at orbits close to 
the ecliptic plane and the information about the latitudinal structure of the SW is hard 
to obtain. There are a few sources of data on the out of ecliptic SW parameters, but only 
one of them from in situ measurements, namely from Ulysses. 

While direct observations of the SW in the ecliptic plane have been collected for 
many years, information on its latitudinal structure had been available only from in- 
direct observations of the cometary ion tails ( Brandt, Harrington, and Roosen||1975 1, 



until radio-astronomy observations of interplanetary scintillation (Hewish, Scott, and 



Wills||1964 Coles and Maagoe||1972[ i and spaceborne measurements of the Lyman-a 



helioglow d Lallement, Bertaux, and Kurt|1985 Bertaux et al.\\99 5S became available. 
To our knowledge these two techniques remain the only source of global, time-resolved 



information on the the solar wind structure. The launch of the Ulysses spacecraft ( Wen- 



zel et al.\l9S9\ improved our understanding of the 3D behavior of the solar wind by 



offering direct in situ observations and a very high resolution in latitude, but a poor 
resolution in timeQ 

The solar wind structure varies in latitude with the solar activity cycle. Knowl- 
edge of its evolution is needed to construct credible models of the heliosphere and 
its boundary regions. With the history of the SW evolution based on a homogeneous 
series of data and retrieved using homogeneous analysis method, one obtains a tool to 
interpret both present heliospheric observations, such as the ongoing measurements of 
Energetic Neutral Atoms (ENAs) by the NASA Interstellar Boundary Explorer (IBEX, 
McComas et a/.|[2009 1 and in situ measurements of the heliospheric environment by 



the Voyagers, and to compare them with the results from past and current long-lived 
experiments such as Solar Wind ANisotropy (SWAN) onboard SOlar and Heliospheric 
Observatory (SOHO) etc. 

We assume that the solar wind expansion is purely radial, its speed does not change 
with solar distance, and density drops down quadratically with distance from the Sun. 
These assumptions are valid for close distances from Sun (r < 10 At/). Their validity 
at farther distances will be considered in the Discussion section. 

In the following text "CR-averaged" data mean Carrington rotation (CR hereafter) 
averaged values, where Carrington rotation period is the synodic period of solar rota- 
tion equal to 27.2753 days ( Franz and Harper|2002| l. Throughout the paper, "adjusted" 
means scaled from the value xq measured at a heliocentric distance ro to the distance 

= 1 AU, where the value specific for r^ is calculated as x^ - xq (ro/rg)^. 



The same latitudes were visited only a few times during the ~ 20-year mission. 
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2. Datasets used 

In our studies we use 3 complementary data sources: in situ in-ecliptic measurements 
from various spacecraft combined in the OMNI-2 collection, in situ out-of-ecliptic 
measurements of solar wind plasma by Ulysses, and remote-sensing radio-observations 
of interplanetary scintillations (IPS), interpreted using the Computer Assisted Tomog- 
raphy (CAT) technique. 

2.1. in situ in-ecliptic measurements collected from various spacecraft 

Solar wind in the ecliptic plane (which differs from the solar equator plane by 7.25°) 
is a mixture of a "genuine" slow solar wind, fast solar wind from coronal holes, solar 
wind plasma from stream-stream interaction regions, and (intermittently) interplanetary 
coronal mass ejections. The features of these components change with solar activity. 

The parameters of the SW vary with the phases of Solar Cycle (SC hereafter) and, 
as recently discovered by [McComas et oL ( 2008[ l, also secularly. The OMNI-2 data 



collection, available at http://omniweb.gsfc.nasa.gov/ (see King and PapitashviU 2005 1, 
gathers SW measurements from early 1960s until present and brings them to a common 
calibration. The absolute calibration of the current version of the OMNI-2 collection is 



based on the absolute calibration of Wind measurements ( Kasper et al. [2006) 1. 



The main source of data is nowadays Advance Composition Explorer (ACE) and 
Wind. Before the Wind era the main datasource was measurements from Interplanetary 
Monitoring Platform-8 (IMP-8), with gaps filled by miscellaneous spacecraft. The data 
from the epoch before IMP-8 are from various experiments and could not be reliably 
brought to a common calibration with the IMP-8/W/nii system because of the lack of 
overlap between the measurement time intervals. Our analysis starts in 1985, when 
IMP-8 was already in operation (see also review by Bzowski et al.\10\2 1 



The solar wind parameters show considerable variations during one solar rotation 
period, with quasi-periodic changes from slow to fast solar wind speed and related 
changes in density. The time scale of the changes of the fast/slow wind streams is 
comparable to the solar rotation period and thus constructing a full and accurate model 
of the solar wind variation as a function of time and heliolongitude is currently not 
feasible because of the lack of sufficient data. 

Here in our analysis, we average-out heliolongitude variations. We start from the 
hourly values of solar wind density and speed available in the OMNI-2 collection, and 
we construct a time series of Carrington rotation-averaged parameters of the solar wind 
with the grid points set precisely at halves of the CR intervals. Small deviations of the 
averaged times from the halves of the rotation periods are linearly interpolated. Thus, 
we develop an equally-spaced time series of the solar wind in-ecliptic densities adjusted 
to 1 AU and speeds, which are presented in Figure [T] and Figure |2] respectively. From 
these series we obtain a time series of solar wind flux calculated as the product of speed 
and density, shown in Figure [3] 

The time interval shown in Figures [T[ |2] [3] starts before the solar maximum in 1990 
and includes the solar minima in 1995 and 2008, as well as the maximum in 2001. 
Neither the density nor speed seems to be correlated with solar activity. The density 
features a secular change ( [McComas et aZ.|2008[ ), which began just before the last solar 
maximum and leveled off shortly before the present minimum. The overall drop in the 
average solar wind density is on the order of 30%. 
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Figure 1. Camngton rotation averages of solar wind density in the ecliptic plane adjusted to 1 AU, calculated 
from the hourly averages from the OMNI-2 database, with vertical lines indicating the three phases of solar 
number density changes during last 20 years that are described in the text. The two horizontal lines mai'k the 
density average values calculated for the intervals 1985 - 1998 and 2002 - present. 
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Figure 2. Carrington rotation averages of solar wind speed in the ecliptic plane, calculated from the hourly 
averages from the OMNl-2 database, with two times intervals marked by vertical lines indicating the short 
time intervals when the drop in speed is registered. They correspond to the minimum phase of solar cycle. 



During the last 20 years we can see 3 phases of the SW density changes: until 1998 
the normal values phase, when the average number density was ~ 8 cm"^ (approxi- 
mately equal to the average value since the beginning of observations), then between 
1998 and 2002 a phase of a rapid decrease that seems to be correlated with the ascending 
phase of solar activity during SC 23, and finally a phase of low density, with an average 
number density reduced to ~ 5.5 cm"^, which lasts until present. These phases are 
marked with vertical lines in Figure [T] The horizontal lines indicate the average values 
before and after the drop. 

The fluctuations of density are generally anticorrelated with speed fluctuations and 
within the low values phase smaller in magnitude than these in the phase before 1998. 
These can be associated with the persistence of coronal holes at equatorial latitudes, as 
convincingly illustrated by|de Tornal(|2011|). 
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Figure 3. Carrington rotation averages of solar wind flux in the ecliptic plane adjusted to 1 AU, calculated 
from the CR-averages of in-ecliptic solar wind speed and density. The vertical lines show the time intervals 
when the slow long-term drop in the flux values occurs. The horizontal lines mark the average values in flux 
described in the text. 



The changes in solar wind flux adjusted to 1 AU are the most pronounced of all 
discussed parameters. The steady and slow decrease that started after 1995 seems to 
continue to the present despite an 8-year plateau between 1999 and 2007 (see Figure[3]|. 
The changes in flux are of the order of 40% from the average calculated from the data 
interval between 1985 and 1995 to the average calculated from the interval 2008.2 - 
2011.5. 

The three phases in the solar wind flux corresponding to the phases pointed out in 
the discussion of the density variations have different time boundaries; they are marked 
with horizontal and vertical lines in Figure[3]in analogy with Figure[T] The speed shows 
multi-timescale variations, but its average value seems to be basically constant in time 
apart from the last 2-3 years, when a small drop of ~ 3% is observed ( jMcComas et al. 



|2008| l. In Figure|2]we indicate only the two time intervals with a small decrease in SW 
speed that are common in time with the lowest part of solar activity in 1997 and 2008. 

2.2. in situ out-of-ecliptic Ulysses measurements 

Ulysses was launched in October 1990 and has been the first and only spacecraft that 
traversed regions of the heliosphere close to solar poles and provided unique samples 
of solar wind. Its orbit was nearly polar, with aphelion of ~ 5.5 AU and perihelion 
~ 1.4 AU, in a plane almost perpendicular to ecliptic and solar equator. The period 
of the orbit was about 6 years, with the so-called fast scan (from south to north pole 
through perihelion), which lasted ~ 1 year and the slow scan (from north to south pole 
through aphelion), which lasted ~ 5 years. 

The heliolatitude track of Ulysses is shown in Figure |4] with superimposed Earth 
heliolatitude position and a solar activity graph represented by the F10.7 cm flux ( |Cov-| 
|ington|[T969[ |Tapping||l987[ |Svalgaard and Hudson||2010| . The heliocentric distance 
of Ulysses for one complete orbit is shown in Figure [5] The spacecraft was launched 
during the solar maximum conditions and during its 20 years life it observed the Sun 
during the whole solar cycle and further until the last prolonged solar minimum. 

The discoveries and findings from the plasma measurements by the Solar Wind 



Observations Over Poles of the Sun (SWOOPS) experiment on Ulysses (Bame et al. 
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Figure 4. Illustration of heliolatitude position of Ulysses (red) and Earth (green) during the time span of 
Ulysses mission (November 1990 - June 2009). The pale blue line is the F10.7 solar radio flux I Covington] 
|1969|[T^pping|1987) , superimposed to correlate variations in solar activity with Ulysses heliolatitude during 
its more than 3 polar orbits. 
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Figure 5. Ulysses heliocentric distance during one of its polar orbits. The solid line corresponds to the fast 
part of the scan, and the dashed line to the slow part of the scan. The horizontal lines mark the 10-degree 
heliolatitude bins used in the analysis. 



1992") can be found in many papers (see e.g. 'Phil lips et a/.||1995 Marsden and Smith 
1997,|McComas et a/.|2000 i. They include the bimodal structure of solar wind ( jMcCo 



mas et aZ. 11998^ 7 which is fast and uniform at mid- and high-heliolatitudes and more 



variable, slower and denser at lower heliolatitudes during low activity of the solar cycle. 
During maximum of solar activity the solar wind is highly variable at all heliolatitudes, 
with flows of slow and fast wind interleaved ( |McComas et a/.|2003| l. 



Detailed studies of the fast solar wind parameters presented by McComas et al. 



(2000 1 revealed latitudinal gradients in proton speed and density of flows from polar 

coronal holes. These gradients are equal to ~ 1 km s 'deg ' and 0.02cm"''deg 

respectively. Ulysses measurements also confirm the existence of the secular changes 
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in S W parameters reported by in-ecliptic spacecraft. [McComas et al. ( 2008[ l showed 
that the global solar wind exhibits a slight reduction in speed (~ 3%), but a much 
greater one in proton density (~ 17%) and dynamic pressure (~ 22%). This result 
was demonstrated simultaneously at low and high latitudes, allowing these authors to 
conclude that the variations were truly global. It was also recorded that the band of 
the slow SW variability extended to a higher latitude during the last Ulysses orbit (see 
McComas et a/.|2006) than during the first one. Ebert et al. ( 2009) reported that above 



t36° heliolatitude the spatial and radial variability in SW parameters remains consistent 
and relatively small, indicating that the fast solar wind plasma flows from the polar 
coronal holes are steady and uniform. 

Ulysses found that the heliolatitude structure of the SW during the two solar minima 
was largely similar (see Figure |6]l, featuring an equatorial enhancement in density with 
the associated reduction in velocity (the slow wind region), and that during solar max- 
imum the slow wind and fast wind from small coronal holes exist at all heliolatitudes 
(see also middle panel in Figure |8]l. However, the region of slow wind seems to reach 
higher in heliolatitude during the last solar minimum than during the minimum of 1995, 
which is much less conspicuous in density (see Figure|6|. 

Here we study the latitudinal variation of solar wind parameters based mostly on the 
fast latitude scans to avoid possible convolution of heliolatitude, time and heliocentric 
distance effects. We leave the data from the slow scans, which covered almost 5 years 
each (i.e. almost a half of solar cycle) for verification of the derived model. Two of the 
3 fast scans (the first and the third one) occurred during the descending phase of solar 
activity and the middle one during the solar maximum conditions (see Figure[4|l. During 
the fast scans Ulysses was close to the Sun (see Figure |5]l, so possible distance related 
effects in the solar wind are small, in contrast to the slow scans, when distance related 
effect may be significant ( Ebert et al. |2009| l. 

The evolution of solar wind speed and adjusted density and flux obtained by Ulysses 
during the fast latitude scans are compiled in Figure [6] where the parameter values are 
averaged over 10-degree bins in heliolatitude. 
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Figure 6. Solar wind speed (upper panel), adjusted density (middle panel), and adjusted flux (lower panel) 
as a function of heliolatitude for the first (blue), second (red) and third (green) Ulysses fast heliolatitude scan. 
The parameters are averaged over the 10-degree heliolatitude bins marked with vertical lines. 
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2.3. Remote-sensing observations of interplanetary scintillations 



Interplanetary scintillation is a phenomenon of producing diffraction patterns on an 
observer's plane by the interferencing radio waves from a remote compact radio-source 
(like quasar) that are scattered by electron density irregularities (fluctuations) in the 
solar wind (see, e.g., Hewish, Scott, and Wills' 1964; 'Coles and Kaufman 1978; Kojima] 
and Kakinuma|1990, Kojima et a/. 1998 , Kojima et aZ. 2007) . The scintillation signal 



is a sum of waves scattered along the line of sight (LOS) to the observed radio-source. 
Most of the scattering occurs at the closest distances to the Sun (the so-called P point, 
see Figure [7]) along the LOS. This is because the absolute magnitude of the electron 
density fluctuations, which are approximately proportional to electron absolute density, 
rapidly decreases with solar distance ( Coles and Maagoe|1972 l. The IPS observations 
are LOS integrated and to provide reliable information on the SW speed they have to 
be deconvolved. 

The magnitude of electron density fluctuations at a given solar distance is correlated 
with the magnitude of local solar wind speed. Thus, inferring the level of electron 
density fluctuations from IPS observations one can estimate the local speed of solar 



wind (Hewi sh, Scott, and WiUs|1964[ [Jackson et a/.|1997[[Kbjima et a/. |1998[ [Jackson 
[eFaL 2003) . However, one needs a formula that links the electron density fluctuations 
6ne with solar wind speed v. Usually, a relation Srie cc v'^ is used. The index y has 
to be established empirically. Such an analysis for various solar wind conditions was 
performed by Asai et aT\\\99^ . 

Observations performed using a system of several radio antennas that are longi- 
tudinally distributed on the Earth allows to infer a very detailed information on the 
structure of solar wind both inside and outside the ecliptic plane ( |Coles and Rickett] 
1976 Tokumaru, Kojima, and Fujiki[|2010| l. To that purpose best suitable seems the 



Computer Assisted Tomography technique (Jackson et al. 1998 Kojima et a/. [[1998 
Hick and Jackson|2001 Jackson et a/. [2003, Kojima et a/. [2007 1. 



The accuracy of the LOS deconvolution depends on (1) geometrical considerations 
(e.g., the geographical latitude of the telescopes, the tilt of the ecliptic to the solar 
equator), (2) the number of observations available, and (3) the fidelity of the correlation 
formula between solar wind electron density fluctuations and speed. While correlating 



these quantities was possible early for the equatorial solar wind ( |Harmon 1975 1, the out 
of ecliptic IPS measurements could only be calibrated once in situ data from Ulysses 



became available (Kojima et al. 2001 1. Such a calibration should be repeated separately 
for each solar cycle because solar wind features secular changes, as discussed earlier in 
this paper. Still, even before the introduction of the CAT technique, the IPS observations 
suggested (e.g., Kojima and Kakinuma|1987 1 that the solar wind structure varies with 
solar activity, with high speeds in the polar regions during low solar activity and the 
slow wind expanding to polar regions when the activity is high. This lends credibility 
to the IPS technique as a source of information on the global heliolatitude structure of 
solar wind speed. 

An extensive program of solar wind IPS observations, initiated in 1980s in the 
Solar- Terrestrial Environmental Laboratory (STEL) at the Nagoya University in Japan 
( Kojima and Kakinuma[[T990| l, resulted in a homogeneous data set that spans almost 
three solar cycles. This data set enables detailed studies of the evolution of solar wind 
speed profile with changes in solar activity ( [Kojima and Kakinuma|198"7{ Kojima et al. 
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Figure 7. Illustration of the geometry of the line of sight in remote-sensing observation of a polar region in 
the solar wind. The observer is close to the Earth in the ecliptic plane and aims its radio-telescope antenna at 
a target so that the LOS crosses a cone with a small opening angle centered at the pole. The signal is collected 
from the full length of the LOS, but the contributions from various parts are different (for details of weighting 
function see |Kojima and Kakinuma|1987^ . In the case of IPS observations, the strongest contribution to the 
signal is from the point nearest to the Sun along the LOS, marked with P, because the source function of the 
scintillation signal drops down with the square of solar distance. Note that the angular area of the polar region 
is quite small, so with an observations program that maps the entire sky only a small region in the map indeed 
includes the signal from polar regions. 
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Tokumaru, Kojima, and Fujik i|20T0 1. 
In the analysis we used data from 1990 to 2011, obtained from the tomography 
processing of observations from 3 antennas (Toyokawa, Fuji and Sugadaira), and in 
addition from another antenna (Kiso) since 1994. The 4-antenna system was operated 
until 2005, when the Toyokawa antenna was closed ([Tokumaru, Kojima, and Fujiki 



|2010| l; since then the system again was operated in a 3-antenna setup. We did not use 
the data collected in 2010, because the number of observations available was too small 
to obtain a reliable result on solar wind structure from the CAT analysis. Since 201 1 the 
STEL IPS observations are again regular with the updated multi-station system. 

The IPS data from STEL are typically collected on a daily basis during ~ 1 1 Car- 
rington rotations per year: there is a break in winter because the antennas get covered 
with snow. Each day, 30-40 LOSs for selected scintillating radio-sources are observed. 
The lines of sight are projected on the source surface at 2.5 solar radii, which is used as 
a reference surface in the time-sequence tomography. 

The latitude coverage of the sky by the time-sequence IPS observations is not uni- 
form and strongly correlated with the Sun position on sky, which changes during the 
year, and with the target distribution on sky. Relatively few of them are located near 
solar poles, because the polar regions are only a small portion of the sky (see Figure|7]l. 
Additionally, the observations of the south pole are of lower quality than these of the 
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north pole because of the low elevation of the Sun during winter in Japan. The original 
latitude coverage was improved owing to the new antenna added to the system in 1994 
and by optimization of the choice of the targets. 

Thus, the accuracy of the remote-sensing measurements of solar wind speed de- 
creases with latitude because of geometry. The polar values are the most uncertain (and 
possibly biased) because the signal in the polar lines of sight is only partly formed in 
the actual polar region of space, which can be understood from the sketch presented in 
Figure |7] 

2.4. Comparison of IPS solar wind speed profiles with Ulysses data 

To verify the results obtained from the IPS CAT analysis, we compared them with the 
data from the three Ulysses fast latitude scans. The Ulysses speed profiles used for 
this comparison were constructed from subsets of hourly averages available from the 
National Space Science Data Center of NASA (NSSDC), split into identical 10-degree 
heliolatitude bins and averaged. They are shown in Figure|8]in red. Since the acquisition 
of the Ulysses profiles took one year each and the first and second scan straddle the 
break of calendar year, we show the IPS results for the years straddling the fast latitude 
scans; they are presented in blue and gray in Figure[8] 

Generally the IPS and Ulysses profiles agree quite well. The sawtooth feature in the 
Ulysses profiles from the second fast scan and in the equatorial part of the first and 
third fast scans is due to the short time Ulysses was sampling the 10-degree bins. The 
fast scans were performed within the perihelion part of the Ulysses elliptical orbit, with 
the perihelion close to the solar equator plane. Hence, the angular speed of its motion 
was highest close to equator and traversing the 10-degree bin took it less than one solar 
rotation period (see Figure |9]l. Thus the sawtooth is an effect of incomplete Carrington 
longitude coverage of the bimodal solar wind by Ulysses, with slow wind interleaved 
with fast wind streams. 

Near the poles the angular speed was slower and it took more than 1 CR to scan 
the 10-degree bin. Thus, when the slow wind engulfed the whole space, the sawtooth 
effect expanded into the full heliolatitude span. By contrast, during the low-activity 
scans the solar wind speed at high latitude was stable, which resulted in the lack of the 
small-scale latitude variations in the CR-averages at high latitudes, despite the uneven 
heliolongitude coverage (see Figure|9]l. 

The IPS yearly averages do not show the short scale latitudinal variability of the 
solar wind speed seen in the Ulysses data (the sawtooth feature) presented in Fig- 
ure |8] because this variability was smoothed by averaging over multiple full Carrington 
rotations, with the alternating fast/slow streams averaged out. 

The difference between the blue and gray lines in the top and middle panels of 
Figure [8] is a measure of true variation of the latitudinal profile of solar wind speed 
during one year Ulysses was going from south to north during the fast latitude scans 
(see Figure |4|i, so the south limb of the profile from Ulysses ought to be closer to the 
southern limb of the blue profile obtained from the IPS analysis, while the north limb of 
the Ulysses profile should agree better with the north limb of the gray IPS profile. And 
indeed, an almost perfect agreement is observed in the second panel, corresponding 
to solar maximum. In our opinion, this is a very interesting observation because it (1) 
shows how rapidly the latitude structure of solar wind varies during the maximum of 



SOLA: solarPhysics_2812-Q4-S4.tex; 5 April 2912; 3:11; p. 11 



Sokot et al. 



SW speed during 1st fast scan 



800 




-80 -60 -40 -20 20 40 60 
heliolatitude [degj 

SW speed during 2nd fast scan 




-80 -60 -40 -20 20 40 60 80 
heliolatitude [degJ 

SW speed during 3rd fast scan 



800 



700 



600 



S. 500 



400 



800 



700 



600 



500 



400 

-40 - 20 20 40 
heliolatitude [degJ 

Figure 8. Solar wind speed profiles from Ulysses measurements and IPS observations for the 3 fast helio- 
latitude scans. Red: Ulysses, blue: IPS during the year of beginning of Ulysses 1st and 2nd fast scan and the 
whole 3rd scan, gray: IPS during the year of the end of Ulysses 1st and 2nd fast scan. Top panel: the first 
fast scan during solar miniimim, iniddle panel: the second fast scan during solar maximuin, bottom panel: the 
third fast scan during minimum. 




solar activity, (2) confirms the credibility of IPS results for solar wind structure, and (3) 
confirms that interplanetary scintillation observations deconvolved by CAT algorithms 
reconstruct the SW structure in a very good agreement with in situ measurements. 

Originally, the interpretation of the speed profiles obtained from Ulysses was not 
clear, it was pondered whether the north-hemisphere increase in solar wind speed was 
a long-standing feature of the solar wind or was it just a time-variability of the wind at 
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Figure 9. Time intervals of scanning the 10 degree bins by Ulysses during fast scans. 



the north pole. Similarly, the question arose whether IPS was able to credibly reproduce 
the solar wind profiles given the fact that some of the profiles obtained approximately 
at the time of the fast scan seemed to disagree with the in situ data (especially the first 
fast scan). The challenging question was which data were more reliable for the yearly 
profile: those from Ulysses that measured features only from the point of solar surface 
it came from, or the yearly profiles from IPS observations that reveal the full-surface 
variations in time? 

The Ulysses data are challenging to interpret in this respect. On one hand the fast 
scans give information on the solar structure from the south to north pole almost in 
one year, but the scans are so fast that the longitudinal variation is convolved with 
the latitudinal structure (see Figure |9]l and on the other hand, the slow scans take so 
long (nearly half of solar cycle) that the parameters obtained for different heholatitude 
bands may differ not only because of a heliolatitude structuring of solar wind, but also 
because of possible time variations of this structure. Without an independent insight 
one is unable to separate between the two. 

Furthermore, the Ulysses heliocentric distances during the fast and slow scans are 
different (see Figure [5]). During the fast scans the distance does not change very much 
(from 2.2 AU close at the poles to 1.4 AU in the ecliptic), so possible distance related 
effects are not pronounced, but during the slow scans, when the distance changes are 
much greater, one can not deconvolve the variation connected with the distance and 
those related to heliolatitude and/or time without a credible MHD model of the solar 
wind behavior between the Sun and the spacecraft. The IPS data provide us with the 
much wanted information on all heliolatitudes simultaneously over one year 

The IPS data are in a very good agreement with OMNI-2 in the ecliptic (see the blue 
and green lines of the corresponding yearly averages in Figure 10 1. Also the agreement 
between the Ulysses measurements obtained from the two Carrington rotations during 
the passages of 0° heliolatitude at aphelion and the OMNI-2 data is almost perfect (cf 
red diamonds in Figure 10 1, which is clearly seen for the passages in 1998 and 2004. 



The apparent discrepancy in 1991 is a result of a fluctuation, as explained in Figure 1 1 



which presents the Ulysses hourly time series compared with the corresponding hourly 
time series from OMNI-2 for the 2 CRs when Ulysses was in the ecliptic plane. The two 



SOLA: solarPhysics_2812-Q4-S4.tex; 5 April 2912; 3:11; p. 13 



Sokot et al. 



650 

600 

550 

c 500 
I 

■B 450 
u 
& 

400 

350 



300 



1990 



Solar wind speed in thie ecliptic plane 



• IPS 

■ OMNI 

♦ Ulysses slow scan 
i Ulysses fast scan 
» Ulysses 10-deg bin 




1995 



2000 
time [y] 



2005 



2010 



Figure 10. Compaiison of yearly solar wind speed values in the ecliptic plane from IPS and OMNI-2 
database with Ulysses fast scan data from the 10-degree bin around solar equator (black triangles; it is the 0° 
bin value from the 10-degree profiles) and Carrington rotation average around the 0° latitude for fast scans 
(red triangles) and 2 Carrington rotation average for slow scans (red diamonds). 
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Figure 11. Comparison of hourly averages of solar wind speed from the OMNI-2 collection (dots) and 
Ulysses measurements (solid line) obtained during two Carrington rotations just after its launch. The inset 
presents a time series of CR-averages of solar wind speed from the OMNI-2 collection to illustrate that the 
average from the Carrington rotation in question happened to be the minimum over a longer time interval. 



time series are in excellent agreement and the big difference between the Ulysses CR- 
average and the yearly-average from the OMNI-2 data exists because the CR-averaged 
speed during this Carrington rotation in question was minimum within the 12 rotations 



included in the yearly average, as illustrated in the inset in Figure 1 1 



The agreement between the Ulysses measurements and the OMNI-2 data set during 
the fast-scan passages is not straightforward to estimate because Ulysses was passing 
the ecliptic latitude bins in a time shorter than a half of the Carrington rotation (see 
Figure [9]) and at heliolongitudes different from Earth's. In this case averaging over a 
full Carrington rotation (shown by red triangles in Figure 10 1 averages out the heliolon- 
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gitude variation of the solar wind, but likely biases the result due to a possible presence 
of heliolatitudinal gradients. On the other hand, averaging over a 10-degree heliolat- 
itudinal bin around the solar equator (black triangles in Figure [T0| mostly eliminates 
the latitudinal gradients, but does not average out the heliolongitudinal modulation. In 
reality, however, the difference between the two averages turned out to be small, similar 
in magnitude to the differences between the yearly averages of the OMNI-2 and IPS 
solar wind speed time series. 

Summing up this section we conclude that the IPS solar wind speed profiles provide 
a reliable insight into the solar wind structure and its evolution with solar activity phase. 
They agree both with OMNI-2 in ecliptic and with Ulysses out of ecliptic for the time 
intervals when they can be compared directly. 

3. Data processing and model description 

Our goal is to retrieve the solar wind speed and density structure and its evolution in 
time and heliolatitude, beginning from the maximum of SC 22, through the minimum 
and maximum of SC 23 (in 1996 and 2001, respectively) until the most recent prolonged 
minimum. We use all relevant and complementary data sets both in and out of the 
ecliptic plane, to construct a homogeneous set of solar wind parameters, which can be 
further used for calculation of ionization losses of ENAs, modeling of heliosphere and 
global interpretation of the Lyman-o' helioglow etc. 

Our procedure takes as baseline absolute calibration of the OMNI-2 data set both in 
speed and density in the ecliptic plane and the absolute calibration of Ulysses measure- 
ments for density and speed out of ecliptic and interplanetary scintillation observations, 
interpreted by the tomography modeling for speed out of ecliptic. Up to now, no global, 
continuous measurements of solar wind density as a function of heliographic latitude 
have been available and this quantity has to be obtained using indirect methods (see 
[Bzowskief a/.|2012| l. 

Tests revealed that an appropriate balance between the latitudinal resolution of the 
coverage and fidelity of the results is obtained at a subdivision of the data into 10- 
degree heliolatitude bins. Concerning the time resolution, ideal would be Carrington 
rotation averages. Regrettably, such a high resolution seems to be hard to achieve 
for two reasons. First, the time coverage in the data from IPS has gaps that typically 
occur during almost four months at the beginning of each year, which would induce 
an artificial 1-year periodicity in the data. Second, the fast latitude scans by Ulysses 
lasted about 12 months and hence differentiating between the time and latitude effects 
in its measurements is challenging. Thus, the full and reliable latitude structure of solar 
wind can currently be obtained only at a time scale of 1 year and this will be the time 
resolution of the model we present. Furthermore, the typical time of solar wind proton 
travel to the termination shock (TS) is about 1 year for ~ 1 keV protons so the adopted 
time resolution is reasonable for the heliosphere modeling. 

3.1. Solar wind speed profiles 

For our analysis we take the solar wind speed available from 1990 to 2011 mapped at 
the source surface on a grid of 1 1 x 360 x 180 records per year, which corresponds to 
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Figure 12. Comparison of the solar wind speed profiles in 2005 obtained directly from IPS (dots) with the 
model with the smoothing function defined in Equation [T] and Equation |3] applied (solid line). The vertical 
lines indicate the boundary bins between different sections of the model function. In 2005 they are: -60, -40, 
10, 40, 70 bins. The numbers above the lines indicate the sections of the model function. 



a series of 11 Carrington rotations. The data are organized in heliolatitude from 89.5° 



North to 89.5° South in the so-called Carrington maps of solar surface (see UMch and 
Boyden'2006Y 



A comparison of the tomography-derived solar wind speed with the in situ measure- 
ments by Ulysses (see Figure |8]l showed that the accuracy of the tomographic results 
depends on the number of IPS observations available for a given rotation. Intervals with 
a small number of data points clearly tend to underestimate the speed. Consequently, 
we removed from the data the Carrington rotations with the total number of points less 
than 30 000. Rotations with small numbers of available observations typically happen 
at the beginning and at the end of the year and at the edges of data gaps. The selection 
of data by the total number of points per rotation constrained the data mainly to the 
summer and autumn months, when all latitudes are fully sampled. 

The selected subset of data was split into years, and within each yearly subset into 19 
heliolatitude bins, equally spaced from -90° to 90°. They cover the second half of Solar 
Cycle 22 and full Solar Cycle 23. We performed the two-step calculation: first Carring- 
ton rotation averaged values of heliolatitudinal profiles and next the yearly averages 
calculated from them. They are shown in Figure 14 It is worth noting that the bin- 
averaged solar wind profiles for specific Carrington rotations have very similar shapes 
to the corresponding yearly profiles with some scatter, which suggests that the general 
latitude structure is stable over a year and changes only on a time scale comparable with 
the scale of solar activity variations. 

Because of the lack of direct data we had to use the speed profiles to infer density 
profiles, as discussed below. Therefore, we decided to construct a function that retrieves 
the SW speed profile for arbitrary time and heliolatitude and smooths the remnant 
variations in the IPS speed profiles, so that they do not bias the inferred densities. An 
additional benefit of such an approach is the ability to determine latitudinal boundaries 
of the fast wind outside the solar maximum phase of solar activity. 
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We smoothed the speed profiles by fitting an approximating piecewise function 



defined as: 



fl3 + ^730 + 
fl4 + b^(p + 
(22 + ^20 + 

as + bs(p + C5 
a\ + bi4>. 



(p34 < ^ < 023 
045 < < 034 

023 < ^ 012 
056 < < 045 
> 012 
< 056 



(1) 



and 



056 < 045 < 034 < 023 < 012 



(2) 



with the additional requirement that the sections of the function connect smoothly, i.e. 
the first derivative of / (0) is continuous at all heliolatitudes. 

This model approximates the SW speed as a function of heliolatitude by linear rela- 
tions in the polar caps and 4 parabolae at midlatitudes and in the equatorial band, which 
all transition smoothly between the neighboring sections. A sketch of the function 



definition is presented in Figure 12 where the splitting into sections is shown on an 
example profile of the IPS data. 

The ordering in Equation [T| is important because those parts of the profile for which 
the number of observations is the highest and which are at the north hemisphere are 
the most reliable. In this way, first of all we fit the function to part 3, that means to the 
north-near ecliptic latitudes, next to part 4, the south-near ecliptic latitudes, next to the 
part for the mid-latitudes and at the end to the polar regions, but with a higher weight 
for the north hemisphere because the number of observations per year is greater in the 
north hemisphere than in the south. 

Under these assumptions the function comes the following form: 



/(0) = 



' A3 + B30 + C302, 

A4 + B4(p + Ci,(fp-, 

A2 + B2<p + C20^ 
A5 + B50 + C50\ 

Ai -hBi0, 
Ag + B(,4>, 



034 < < 023 
045 < < 034 
023 < < 012 
056 < < 
> 012 
^56 



'45 



(3) 



< I 



where the formulae for the coefficients A,, B, and C, are defined in Table[T] The number 
of free parameters is reduced from 16 in Equation[T]to the following 6: 03, b4,C4, bi,b^, C3.I 
The boundaries 0i2, 023, 034, 045, 056 separate the heliolatitude pieces of the fitted func- 



tion (see Figure 12 1 and are chosen separately for each of the yearly profiles so that the 



fit residuals are minimized. In Table |2] we present the values of the parameters of the 
smooth function for each yearly profile. The residuals of the fits are typically 4% and 
do not exceed 10%. 

Table [3] presents the heliolatitude boundaries adopted for each yearly profile of the 
IPS SW speed. They are also shown in Figure 13 which illustrates how the heliographic 
latitude bands change with time and solar cycle phases and how the ranges of solar wind 
regimes evolve. During the interval 1998-2006 the low-latitude bands of the slow and 
variable SW expand to midlatitudes. The behavior of the boundaries fitted for the south 
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Table 2. Values of the free parameters of the function defined in Equation [T] and Table [T] for all 
years. 



year 


^3 




C4 




f 6 


^3 


1990 


418.566 


1.45736 


0.0337921 


-2.17352 


-3.51384 


-0.0392348 


1991 


483.579 


-0.66172 


0.00308865 


2.55869 


-2.23751 


0.0941952 


1992 


449.381 


-2.06557 


0.138359 


0.614671 


-1.89668 


0.121885 


1993 


429.457 


-1.69448 


0.333417 


1.60133 


-0.0603855 


0.0151053 


1994 


506.33 


-3.77307 


0.00505828 


0.0313874 


-0.93796 


0.319795 


1995 


462.637 


-4.29315 


0.137468 


0.573251 


-0.17428 


0.505656 


1996 


347.645 


1.71087 


0.662762 


0.717709 


-0.719594 


-0.246573 


1997 


348.65 


0.0918378 


0.523877 


0.636597 


-0.6109 


-0.115214 


1998 


446.188 


-0.0377569 


0.0900351 


0.0607183 


-0.106951 


0.112234 


1999 


426.63 


-0.115668 


0.0423373 


2.58806 


-2.32994 


0.0152633 


2000 


452.346 


-0.158137 


0.00223068 


-1.74742 


1.23412 


0.0136059 


2001 


452.326 


-1.07672 


-0.00476791 


2.47585 


-2.29808 


0.113553 


2002 


459.426 


-0.647343 


-0.000129429 


2.0322 


0.0351119 


0.0695614 


2003 


529.919 


-2.50005 


-0.0393349 


0.927344 


-1.94706 


0.0850843 


2004 


453.411 


0.494147 


0.0831141 


1.80231 


-2.96517 


0.058379 


2005 


480.172 


0.705444 


0.0966926 


0.925081 


-0.803117 


0.0732095 


2006 


421.171 


1.1226 


0.217231 


-0.0220226 


-1.3036 


0.0698785 


2007 


480.365 


-3.25418 


0.112857 


-0.811889 


-0.302357 


0.222071 


2008 


519.668 


-2.13859 


0.0872578 


1.39421 


0.255003 


0.452983 


2009 


391.003 


0.523058 


0.336663 


-1.08514 


-1.39841 


-0.0276554 


2011 


477.692 


0.704998 


0.0765034 


3.32028 


0.267478 


-0.0960392 



hemisphere is somewhat different, which is probably due to a lower quality of the IPS 
speed profiles in the south hemisphere because of the observations conditions discussed 
above. 

This model function was used to smooth the yearly profiles obtained from the IPS 
tomography data. We use them in further analysis keeping the 10-degree resolution in 
latitude. As it is seen in Figure 14 the smoothing procedure works very well for all 
years, with slightly higher residuals at higher latitudes. The model can be used equally 
well for the solar minimum and maximum conditions and it can be applied to both IPS 
and Ulysses data (see Figure 15 1. 
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Table 3. Boundaries between sections of the 
smoothing function defined in Equation[T] 



year 


056 


045 


034 


023 


012 


1990 


-70 


-40 


10 


40 


70 


1991 


-70 


-40 


10 


40 


70 


1992 


-50 


-20 


10 


40 


70 
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1994 
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-20 





20 


60 


1995 


-40 


-20 





20 


40 


1996 
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40 


70 


1997 


-60 


-10 
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40 


70 


1998 


-60 


-40 


10 
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70 


1999 
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-50 





50 


70 


2000 
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-50 
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70 


2001 
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2003 
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-30 
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50 


70 
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Figure 13. Boundaries between the sections of the model function defined in Equation [T] and Equation |3] 
fitted to the yearly solar wind heliolatitude profiles obtained from the IPS observations. See Table [5] for 
numerical values. 
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Figure 14. Heliolatitude yearly profiles of solar wind speed. The dots represent the IPS yearly data obtained 
from the CAT analysis fTokum aru, Kojima, and FujikipOlO) and the line shows the smoothing function 
defined in Equation [3| for each year separately. The red points and lines superimposed on the profiles (in 
analogy with FigurelSj are parts of the Ulysses fast scan latitude profiles for years 1994, 1995, 2000, 2001 
and 2007 averaged over 10-degree heliolatitude bins. The average times for the profiles are indicated in the 
panel headers. 
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Figure 15. Results of application of the smootliing function defined in Equation[3]to solar wind speed data 
obtained from the Ulysses fast scans, averaged over 10-degree heliolatitudes bins. The boundaiies between 
the model function sections are: (-50, -30, 0, 20, 40) for the first fast scan, (-70, -20, 0, 50, 70) for the second 
slow scan and (-60, -30, 10, 30, 60) for the third fast Ulysses scan. Blue dots and lines: Ulysses original data, 
gray solid lines: results of the smoothing. 



The linear behavior of the fast solar wind as a function of heliolatitude was dis- 



covered by [McComas, Elliot, and von Steiger ( 2002 1 based on in situ measurements 



by Ulysses. In Figure 16 we show the behavior of SW speed for polar regions for the 
first and third Ulysses polar orbit. There is a slight drop in speed at the north hemi- 
sphere between the first and third scan, while the drop in density is well seen for both 
hemispheres. 

3.2. Solar wind density profiles from density-speed correlation 

There have been no direct solar wind density measurements apart from the Ulysses 
in situ data and in lack of remote-sensing observations (e.g. density deconvolved from 
remote-sensing observations of Lyman-a helioglow from SWAN/SOHO) the density 
has to be estimated using indirect approximate methods. The SW density adjusted to 1 
AU can be approximately inferred from the speed profiles, but the correlation between 
speed and density for different heliolatitudes is challenging to obtain. Some insight was 
provided by Ebert et al. ( 2009]), but in their approach the data from the fast and slow 



scans were treated collectively and it was hard to deconvolve the temporal and spatial 
eff^ects (see Figure |5] for the correlation between the solar distance and heliolatitude). 

Here we decided to use an interim solution and retrieve solar wind density from 
correlations between the solar wind speed and density obtained specially for this project 
from the Ulysses fast latitude scans. During solar minimum the solar wind speed and 
density are correlated in heliolatitude but the correlations seem slightly different be- 
tween the first and third latitude scans because of the secular changes observed in the 
solar wind ( |McComas et a/.|2008| l. The correlations for the fast scans are illustrated in 



Figure 17 We assume the following linear relation between the SW proton density «p 
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Figure 16. Solar wind speed (upper row) and density (lower row) variations near the polar regions during the 
solar minimum conditions from Ulysses measurements. Blue lines represent the beginning and end of the first 
fast scan, green lines the beginning and end of the third fast scan. Solid lines are for Ulysses approaching to 
the pole and the dashed lines are for Ulysses receding from the pole. The dots represent the 1-degree latitude 
averages. 



and speed Vpi 

^p, Ulysses (^p) ~ ^scan ^ scan (4) 

where fl.5ca«> b^can fit separately for the speed and density values averaged over 10- 
degree bins using the ordinary least squares bisector method (see OLS bisector in |Isobe| 
et aZ.|1990^ , which allows for uncertainty in both ordinate and abscissa, separately for 



the first and third latitude scans. For the first scan (blue line and points in Figure 17 1 we 
obtain: 

Qfirst = 12.69 + 1.17, bfirs, = -0.01332 + 0.00154, (5) 

and for the third scan (green line and dots in Figure [TTji the correlation formula param- 
eters are: 

a,hird = 10.01 ± 0.65, b,hird = -0.01 107 ± 0.00094. (6) 

Thus, the slopes are almost identical and the main difference between the two formulae 
is in the intercept, which reflects the overall secular decrease in solar wind density 
between the two solar minima. 

The relation between density and speed for the second scan, which occurred during 
solar maximum, is hard to establish directly because there are few points with high 
speed values. In this case the spatial and temporal effects seem to be convolved (as 
discussed earlier). Therefore, we propose to use a very simple relation based on the 
assumption that the density-speed coiTelation does not change with solar activity phase 
and use an arithmetic mean of the relations for the first and third scans as an valid for 
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Figure 17. Correlation between the solar wind adjusted density and speed obtained from the Ulysses fast 
latitude scans. Blue corresponds to the first scan (see Figure |4j, green to the third scan and red to the 
second scan, performed during the solar maximum conditions. The dots correspond to the density-speed pairs 
averaged over the 10-degree heliolatitude bins, the blue and green lines are the linear coirelations specified in 
Equation|4] The dotted red line is the density-speed relation proposed for the transition interval close to the 
solar maximum in 2002, calculated as an arithmetic mean of the correlation relations obtained from the first 
and third latitude scans. 
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Figure 18. Comparison of adjusted solar wind density from in situ measurement by Ulysses during the 
second fast scan during solar maximum (red) with the density calculated from the Equation[4]and Equation^ 
adopted to speed values from Ulysses data from second fast scan (gray). 



the middle time period: 

asecond = [a first + aihird) /2 and bsecond - {b first + btlvrd) /2. (7) 



Relation from Equation |7] is shown in Figure 17 as the red broken line. It seems to 
reconstruct the SW density reasonably well, as shown in Figure 18 where a compar- 
ison of the density values actually measured during the second fast latitude scan and 
calculated from the adopted correlation formula is presented. 
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Since we have three different correlation formulae, we have to specify the time in- 
tervals of the applicability of them. We connect them with the changes in SW density in 
time and adopt the following rules: the formula from the first scan applies to the interval 
before 1998, because before this time no long-term density changes were observed; the 
second relation to the interval from 1998 until 2002, when the density decrease was the 
most visible; and the relation from the third scan for the interval since 2002, when the 
density seems to have changed its temporal gradient. 

The correlation formulae mentioned seem to have a purely statistical character and 
are only good to reproduce the large-scale relation between bin-averaged speed and 
density. They do not allow to reliably reproduce density short-scale variations within 
the equatorial solar wind because, as we verified, there is no significant correlation 
between long-term averages of density and speed in this regime. 

We calculate the yearly profiles of solar wind density as a function of heliolatitude by 
applying Equation |4] to the smoothed speed profiles reported in the preceding section. 



The results are shown in Figure 19 Despite the limitations of the correlation formulae 



the agreement with Ulysses profiles is quite good, especially for polar regions during 



the third fast scan and the northern limb during the first fast scan (see Figure 22 1. There 
is a slight difference between Ulysses measurements and model density results for solar 
maximum in the south polar region in 2000 and 2001, which might be due the following 
reasons: the profile presents the yearly average value and Ulysses the current conditions 
on the Sun; the data are for solar maximum, when the conditions on the Sun change 
dynamically in short time scales; the yearly value from IPS is calculated only from 
about ~ 8 CRs without winter months and Ulysses sampled the highest south latitudes 
just in winter months; and the correlation formula adopted for the solar maximum is 
retrieved from the solar minimum conditions and thus may be not fully adequate. 

3.3. Interpolation in time 

The last step to retrieve the solar wind structure (both in speed and density) as a function 
of heliolatitude and time is a Unear interpolation to nodes at halves of Carrington rota- 
tions. To be consistent with the direct measurement in the ecliptic plane we replace the 
equatorial bin obtained from the presented analysis with the CR-averaged time series 
from OMNI-2, linearly interpolated to halves of CRs. The ±10° bins are replaced with 
values linearly interpolated between the ±20° bins (respectively) and the equatorial bin. 
Because the data from Ulysses are available only to ~ 80° heliolatitude and the IPS data 
at ±90° bins are scarce and thus not fully reliable, we calculate the polar values from a 
parabolic interpolation between the ±70° and ±80° bins. We fit a parabolic curve to the 
points from bins ±70° and ±80° and their mirror reflection around the appropriate pole 
and we calculate the values for the ±90° heliolatitude from the fits. 

As a result of such a treatment, we have utilized all available information on the 
equatorial band of solar wind traversed by the Earth during its yearly orbital motion, 
and away from the equatorial band, where such detailed information is not available, 
we have a smooth transition into the region of the low time-resolution model. At the 
poles, because of the extra smoothing procedure, we avoid a conical sharp peak, having 
latitudinal gradient precisely zero. 

The one-year gap for 2010, when data on solar wind speed from IPS analysis are not 
sufficient to retrieve solar wind speed, we fill in by the average value calculated from 
the straddling years 2009 and 201 1. 
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Figure 19. Heliolatitude yearly profiles of solar wind density obtained from the density-speed coirelation 
discussed in the text adopted to the smooth speed profiles obtained in the previous subsection. The red dots 
and lines show the corresponding parts of the Ulysses fast scan profiles, similarly as in Figure[T4] The average 
times for the profiles are indicated in the panel headers. 
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4. Results 



The primary results of this study are presented in Figure 20 which shows heliolatitude 
vs time maps of solar wind speed and number density adjusted to 1 AU. The results 
confirm that solar wind speed is bimodal during solar minimum, slow at latitudes close 
to solar equator (and thus to the ecliptic plane) and fast at the poles. 

The heliolatitude structure evolves with the solar activity cycle and becomes flatter 
when the activity is increasing. The structure is approximately homogeneous in helio- 
graphic latitude only during a short time interval during the peak of solar activity, when 



solar wind at all heliolatitudes is slow (see the panel for 2000 in Figure 14 1 and highly 
variable. 

Shortly after the activity maximum the bimodal structure reappears and the fast wind 
at the poles is observed again, but switchovers from the slow to fast wind close to the 
poles may still occur during the high activity period: compare the panels for 2001, 2002 
and 2003 in the aforementioned figure and see the solar activity level depicted with the 
blue line in Figure |4] 

During the descending and ascending phases of solar activity there is a wide band 
of slow and variable solar wind that is on both sides of the equator and extends to 
midlatitudes. The fast wind is restricted to polar caps and upper midlatitudes. 

At solar minimum, the structure is sharp and stable during a few years straddling 
the turn of solar cycles, with high speed at the poles and at midlatitudes and a rapid 
decrease at the equatorial band. Thus, apart from short time intervals at the maximum 
of solar activity, the solar wind structure close to the poles is almost flat, with a steady 
fast speed value typical for wide polar coronal holes, which is in perfect agreement with 



the measurements from Ulysses ( Philhps et a/.||1995| McComas, Gosling, and Skoug 
[20001 IMcComas et a/.|2006| l. 

The bimodal structure of solar wind speed is also reconstructed in solar wind density. 
The variable dense flows are at low latitudes and rarefied near the poles. The solar wind 
density changes are anticorrelated with speed. The dense plasma flows are recorded at 
all heliographic latitudes only during the peak of solar maximum phase (see 2000 in 



Figure 19 1. For other years the solar wind with a low number density appears at higher 
latitudes with the minimum values during solar minimum. 

Figure [19] also shows that during minimum of SC 23 the structure of density was 
more narrow around equator than it was during the last minimum of SC 24. It seems 
that the slow and dense plasma flows typical for solar minimum conditions extend to 
higher latitudes (about 10° farther) than it was during the previous cycle. It means that 
the secular changes in solar wind density are very well reflected in our results (the 
profile width at midlatitudes is wider during the most recent years). 

The maps of solar wind speed and density in Figure |20] show also a slight hemi- 
spheric asymmetry, that seems to reverse from one SC to another 

The evolution of solar wind parameters as a function of time and heliolatitude is 
needed in the modeling of the heliosphere and in the interpretation on measurements 
such as ENAs observations by IBEX (M cComas et a/. [20 09 ) and heliospheric Lyman-a 
glow analysis ( jQuemerais et a/.|2006[|Lallement et a/.|2010a| l. Having the density and 
speed we can easily calculate solar wind flux: 

F [(Pj, f,) = Vp [cf)j, t,) Hp [(f>j, ti) , (8) 
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and charge-exchange rate between solar wind protons and stationary H atoms: 

Pcx [4>}, t) = Hp (0;, t) Vp f,) crcx (vp r,)) , (9) 

where crcx is the cross section for charge-exchange rate ( [Lindsay and Stebbings|2003] ), 
and dynamic pressure: 



(10) 



where j numbers the 10-degree heholatitude bins and / the number of Carrington rota- 
tion in our time grid. 



Figure 21 shows contour maps of solar wind flux, charge-exchange rate and dynamic 
pressure for the years since 1990 to the end of 201 1 and for heholatitude from -90° to 
90°. The gap in 2010 is filled by the average value calculated from 2009 and 201 1. The 
flux features a clear secular drop after the last solar maximum. The bimodal structure 
during the current solar minimum seems to be even better defined than during the pre- 
vious one. The structure at solar maximum is quite flat and seen longer than in the case 
of speed. 
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The charge exchange rate basically follows the behavior of the flux, with a clear lat- 
itudinal contrast during low activity period, an almost flat structure at solar maximum, 
and the secular drop after the most recent solar maximum. 

Dynamic pressure behaves different than the flux. While variations in time are clearly 
visible, the latitudinal structure is much less pronounced and does not vary much during 
the solar cycle. Basically, dynamic pressure is almost spherically symmetric (with pos- 
sible exceptions in the polar caps, which, however, cease closer to the poles than in the 
case of the flux), and the most striking feature is the secular drop in the strength, which 
begins earlier than in the flux, namely about 1998. Such a behavior has pronounced 
consequences for the shape of the termination shock, which should not feature a very 
strong latitudinal variation in size, but which should now be significantly closer to the 
Sun than during the previous solar minimum. 



5. Discussion 

To check the credibility of our results we compare them with the Ulysses data from all 
scans. The Ulysses data are prepared by splitting the hourly data into full Carrington ro- 
tations and calculating average values. Next we linearly interpolate our model values to 
the times and heliolatitudes corresponding to the Carrington rotation-averaged Ulysses 
data. 

A comparison is presented in Figure |22] The agreement is quite good in the ecliptic 
parts of the Ulysses orbit and satisfactory for higher latitudes. The model retrieves 
the fast solar wind speed, but some discrepancies exist for the slow and variable solar 
plasma flows. We have a better agreement in solar wind speed than in density, which 
is understandable, since the density values were derived from the already approximate 
speed values using approximate statistically-derived solar wind density-speed relations. 

The worst agreement between Ulysses density measurements and our model is dur- 
ing the third slow scan (during descending phase of solar activity), when Ulysses in situ 
measurements for this solar cycle phase give highly variable values, nearly 50% greater 
than our model predictions. The source of the disagreement might be connected with the 
density-speed correlation formula we adopted for this time interval: for solar maximum 
we use an average formula from the two fast scans during solar minimum assuming that 
the correlation does not change with the solar cycle. 

The overall agreement, however, is much better The residuals of speed are typically 
10%, not exceeding 30%, and typical residuals of density are 20 - 30%, not exceeding 
60%. The sign of the residuals varies, which suggests there is no systematic global bias 
in our method. Given all the uncertainties in the absolute calibration of both in situ 
measurements and IPS observations and the relative simplicity of our approach, we 
believe such a level of agreement between the model and the measurements is quite 
good and hard to improve on without an additional data source. 

Such a source of additional information might result from an inversion of pho- 
tometric maps of the Lyman-or helioglow obtained from SWAN/SOHO observations 
( |Lallement et a/.|[2010b[ ), aimed at calculation of the total ionization rate of neutral 
interstellar hydrogen in the inner heliosphere as a function of heliolatitude and Carring- 
ton rotation. With this information, one might be able to follow the idea presented by 



Bzowski et al. ( 2012| l and independently calculate the profiles of solar wind density. 
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Figure 21. Contour maps of the solar wind flux (in 10* cm ^ s '), charge-exchange (in 10 ' s ') and 
dynamic pressure (in nPa) shown as a function of time and heliolatitude. 



In this analysis we assumed that the solar wind parameters obtained from OMNI-2, 
Ulysses and IPS are directly comparable, i.e. that there is no systematic change in solar 
wind speed with the solar distance between the region from the solar wind acceleration 
region (a dozen solar radii) to 1 AU, from which the IPS solar wind speeds are retrieved, 
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and Ulysses, which measured between ~ 1.4 and ~ 5.5 AU. Also, we assume there is 
no distance-related change in the density other than the simple 1/r^ scaling. 

The assumption of purely radial expansion of solar wind implies that the latitude 
structure we inferred in this paper does not change until the termination shock. This 
may not be exactly true, as pointed out by ( |Fahr and Scherer 2004 1, who argue that 
the pickup ions (PUIs) pressure induces nonradial flows at large heliocentric distances. 
Such flows would cause changes both in the radial component of solar wind speed 
and in the local density. A more thorough study of this eff'ect requires, in our opin- 
ion, using a multifluid, 3D and time dependent MHD modeling with turbulence and 
boundary/initial conditions taken from observations. Such a model, to our knowledge, 
is still in development ( Usmanov et a/.|20TT i. Our results seem to be well suited as the 
boundary conditions for such modeling. 

The assumption of a drop in density and constancy of speed is increasingly 
invalid with an increase of solar distance because of the interaction with neutral inter- 
stellar gas, which results in the creation of pickup ions and slowdown and heating of the 



distant solar wind. These phenomena were extensively discussed by Fahr and Rucinski 



(fT999ll; iFdir and Ruciriski] pOOT] [2002) 1; |Fair| ( |2007] l; |Lee et al.\ ( |2009| l; IRicharisoii 
[eFaT] ( fT995ir2008b|a| l. 

However, for the global modeling of the heliosphere and calculation of survival 
probabilities it is the total flux of solar wind, being a sum of the core solar wind and 
pickup ions, that counts most. In this respect, the total solar wind flu x quite exactly 
follows the scaling due to the continuity conditions, as discussed by Bzowski et al. 



( 2012| l. This is because a vast majority of PUIs are created due to charge exchange and 
thus are not new members of the solar wind population. 
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Figure 22. Comparison of Ulysses in situ measurements for slow and fast scans with model results for solar 
wind speed and adjusted density. Presented are Carrington rotation averaged data. Blue - model, red - Ulysses. 
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6. Summary and outlook 



We have combined in situ and IPS remote-sensing observations to retrieve the heholat- 
itude structure at 1 AU of solar wind speed and density and its evolution in time and 
produced a time series of heliolatitude profiles of solar wind speed and density for each 
year between 1990 and 2009, averaged over 10-degree heliolatitude bins. 

We carefully checked the agreement between the solar wind speed measurements 
available in the OMNI-2 data base and Ulysses fast latitude scans and IPS measure- 
ments of solar wind speed in the ecliptic plane and found it to be excellent. We verified 
the agreement between the yearly-averaged heliolatitude profiles of solar wind speed 
available from Computer Assisted Tomography processed the IPS observations and 
the speed profiles obtained from the three fast latitude scans from Ulysses. Thus, we 
adopted the IPS -derived yearly speed profiles with some additional smoothing as rep- 
resentative for solar wind starting from the maximum of solar activity in 1990 until the 
end of available data in 201 1 (see Figure 15 i. 

Based on the measurements from Ulysses obtained from the first and third fast lati- 
tude scans, performed during the previous and most recent solar minima, we established 
approximate linear correlations between solar wind density and speed that can be used 
only to retrieve the heliolatitude structure of solar wind density. We found that the slope 
was practically unchanged between the two cycles, but the intercept changed because 
of the global reduction in solar wind density observed since the solar minimum in 2001. 
Using these correlations, we calculated the yearly solar wind density profiles based on 
the smoothed velocity profiles (see Figure [T9|. 

To facilitate the use of our time and heliolatitude model of solar wind structure in the 
interpretation of heliospheric measurements, we calculated bilinearly interpolated pro- 
files of solar wind speed and density on a Carrington rotation period grid and replaced 
the equatorial values obtained from the aforementioned procedure with Carrington ro- 
tation averages of the solar wind speed and density available from the OMNI-2 time 
series (see Figure 20 1, and +10° and ±90° bins by appropriate interpolation values. 

Further, we calculated some quantities frequently used in the heliospheric studies, 
including solar wind flux, charge-exchange rate with neutral H, and solar wind dynamic 



pressure (see Figure 21 



The results presented in this paper under the assumption about the radial solar wind 
flow can be applied in global heliospheric modeling, where on one hand one needs to 
know the long-term variations in solar wind ram pressure and production rate of pickup 
ions, but on the other hand short-scale variations in the parameters are less important. 
They can also be applied in the interpretation of in-ecliptic heliospheric measurements 
such as observations of H ENA flux by IBEX or observations of the Lyman-c helioglow, 
where most important are effects of the solar wind interaction with hydrogen atoms 
within ~ 10 AU from the Sun. A homogeneous treatment of long time series of solar 
wind observations enables a direct application of our results to the interpretation of 
heliospheric experiments from 1990 until present. 

The model can still be improved by additional sources of data and by a more ad- 
vanced modeling of solar wind evolution in time and solar distance. Improving the 
model will most likely have to be an iterative process. For example, the current model 
could be used as input to a procedure used to fit a ionization rate model to the global 
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helioglow observations by SWAN/SOHO, and from the result of such fitting an approx- 
imation of solar wind density profiles and their evolution in time could be retrieved and 
used to improve our solar wind evolution model. Studies of the time and heUographic 
latitude behavior of solar wind is thus still a work in progress. 
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